% For "The Climate Risk Premium" (Lemoine, JAERE, 2020)

% sub_damagenodes:
% Hard-code quadrature nodes for the damage distribution, as calculated from Fortran90 version of truncated_normal_rule.

switch option_damagetype
    case 'base'
        switch Params.quadrature_nodes
            case 8
                if option_ocean~=1
                    nodes_damage =   [        -13.68597003642037
                        -11.86271564331279
                        -10.32361097414228
                        -8.933801551182290
                        -7.655041961042801
                        -6.494008164928736
                        -5.509542684524233
                        -4.841957685025910       ];
                    weights_damage = [       0.8351137118378431E-05
                        0.1088826414072052E-02
                        0.2162329874334853E-01
                        0.1280148943376288
                        0.2969047990109475
                        0.3169173815542833
                        0.1787272281973424
                        0.5671522060525946E-01       ];
                else
                    nodes_damage =     [ -13.51405607266014
                        -11.69172422636591
                        -10.15334923061702
                        -8.764119272668747
                        -7.485762302816036
                        -6.324882439454035
                        -5.340238761293922
                        -4.672232347190501  ];
                    weights_damage = [   0.8286502766903319E-05
                        0.1082482800438573E-02
                        0.2153607165348518E-01
                        0.1277205212343216
                        0.2967047195513303
                        0.3171350610752351
                        0.1790001687927157
                        0.5681268838970623E-01 ];
                end
            case 9
                if option_ocean~=1
                    nodes_damage =   [        -14.25519027482930
                        -12.47705872602514
                        -10.98080907389941
                        -9.631855947663674
                        -8.387182411113939
                        -7.239095804416475
                        -6.209552836202054
                        -5.360912782346566
                        -4.808465149412803       ];
                    weights_damage = [       0.1207931844578181E-05
                        0.2252984476218343E-03
                        0.6361621684156336E-02
                        0.5446907893735083E-01
                        0.1875694256549366
                        0.3047509602740724
                        0.2641109043172460
                        0.1380063908640963
                        0.4450511188867488E-01     ];
                else
                    error('Not hard-coded')
                end
            otherwise
                error('Not hard-coded')
        end
        
    case 'level'
        
        % 8 quadrature nodes
        nodes_damage =   [      -10.18805673710705
            -8.586977549712746
            -7.232763642389899
            -6.006531755449762
            -4.873594187168683
            -3.838230744532030
            -2.951247267481227
            -2.341237073455048      ];
        weights_damage = [          0.9858324624714633E-05
            0.1255877191644031E-02
            0.2428517162222100E-01
            0.1391218490931034
            0.3093168937951603
            0.3124835081440694
            0.1646032861830343
            0.4892355564614301E-01        ];
end

